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O ■ Abstract 
00 

Q> ■ We study the relativistic van Alphen-de Haas effect in cold dense quark matter in 

r~| ' a homogeneous magnetic field which arises from oscillations in the effective potential 

Q-ii of an underlying extended NJL-model. Finally, we discuss the phase structure of the 

Q^' model in the parameter plane (/i, H) and the solution of the mass gap equation. 

1 Introduction 

X. 

^ , The study of dense matter phenomena is an important branch of physics. Most properties 
of matter are only weakly affected by changes in external parameters, such as temperature, 
chemical potential, electromagnetic field etc., when they vary within ranges far from extreme 
values, i.e. under laboratory conditions. Therefore, each time when a small change in 
external parameters leads to a considerable response is remarkable and means a discovery of 
a significant effect. We can mention as examples the quantum Hall effect and oscillations of 
the magnetic moment predicted in 1930 by Landau and experimentally observed by van 
Alphen - de Haas 0. Sixty years have passed since the discovery of this effect, but untill 
this day only quantum oscillations considered in the nonrelativistic case have been discussed 
extensively. Most attention of researchers dealing with oscillations of the magnetic moment 
is now focused on the relativistic problem since the results of these studies may be applied 
to cosmology, astrophysics and high-energy physics P, 
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The consideration of dense quark matter as described by a chemical potential leads to 
screening of quark charges and possibly to a restoration of chiral symmetry. In order to 
describe the corresponding chiral phase transition we must have a realistic model. It is 
well known that models of the Nambu-Jona-Lasinio (NJL) type with inclusion of back- 
ground fields, temperature and chemical potential are a good laboratory for investigating 
the nonperturbative phenomenon of spontaneous breaking of chiral symmetry, as well as for 
describing the low-energy sector of quantum chromodynamics p, |^, 0] 

In this paper we shall in particular demonstrate the possibility of the existence of the 
van Alphen - de Haas effect in quark matter with finite density in a homogeneous constant 
magnetic field. As model Lagrangian for light quarks we shall use the NJL model where 
quarks acquire a dynamical mass due to spontaneous breaking of chiral symmetry in the 
presence of a homogeneous magnetic background field. 



2 Effective potential of the NJL model 

The Lagrangian of the NJL model with light quarks u, d of the flavor group SU (2) / and the 
color group SU (Nc) has the form (flavor and color indices of quarks will be suppressed) 

L = qil^D^q + ^[{qqf + {qi^'>Tq)\ (1) 

where the covariant derivative is given by = — iQA^, Q = diag{ei, 62), G is a four- 
fermion coupling constant, and r are isospin Pauli matrices. As is well known, this model 
is symmetric under continuous chiral transformations. 

To investigate the properties of the vacuum of the NJL model it is convenient to employ 
instead of the quark Lagrangian (|Tp the quark-meson Lagrangian 

= qi^^D^q-qia + t^'T7f)q-^J:^ (2) 

where = o"^ + vf^, which, due to the equations of motion for the fields a, tt, is equivalent 
to (|l]). It is convenient to consider the effective action of the model which is given by 

exp{zS'e//(/i, 0-, tt)} = J DqDq exp (^i J {L„ + yUg7og) d'^x^ , 

where 

^S^ff = - I d'x^^^^^^ - iTr HtYD, + /X70 - a - t^'fn), (3) 

and fi denotes the chemical potential. 

Let us first recall the well-known properties of the model with vanishing chemical po- 
tential and vanishing background magnetic field. Assuming, as usual, that the fields a, vf do 
not depend on the space-time, and using the definition 
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we obtain the effective potential 



^^e// = ^ + 4zy^ln(S -P) = ^. (4) 

Switching in V"e//(cr, tt) to the Euchdean metric and introducing a cut-off for the integration 
region < A^, we obtain 

^\/o(S) = SV2G'-(87r2)-2{A4in(l + sVA2) + 

+ A2s2-EMn(l + AVS')}. (5) 

The condition that the function (|^) is stationary is 

(9Vo(S)/9S = 0. (6) 

Obviously, for G < Gc = 47r^/A^ equation (|]) has no solutions except the trivial one S = 0. 
On the other hand, for G > Gc a. nontrivial solution S = (Tq does exist (we choose ttq = 0, 
which means that a pion condensate is absent), and the potential Vo(S = cxo) has a global 
minimum at the point S 7^ 0, signalling spontaneous breakdown of chiral symmetry and the 
appearance of a quark mass. 

In the following, we shall study a gas of quarks with non-zero chemical potential in 
the presence of an external homogeneous magnetic field which is determined by the vector 
potential = H5^2Xi. In this case the energy spectrum of quarks is given by 

El^ = T?+pl + 2\e,\Hn, (7) 

where n = 0,1, 2,... is the Landau quantum number. Let us next temporarily assume that, 
together with the chemical potential fi, a heatbath with temperature T acts on the system 
described by the NJL Lagrangian In this case, to obtain the effective potential the 
following transformation of the integration variable must be made — ^ 27r(/-|-l/2) /f], I = 
0, ±1, ±2, where 1/(3 = T. After that the effective potential of our model can be written 
as 

where V^/3(S)/A^c denotes the part of the effective potential with 7^ and H = 0, and 
Vo/3(S, \ei\H)/Nc is the part with ^ = 0, H ^ 0. Finally, AV{fi,P, S, \ei\H) is responsible for 
the situation in which both fi and H are non-zero. The third part of the effective potential 
is of particular interest here and takes the form 



-^AV^(/i,/?,E,|e.|i/) = -E^^/ dpsY.(^n{Hi^ + e-^'-^^^-+^^)x 

i=i vr^ Jo 

(1 + e-/3(s,.„-M))]}^ (9) 
a„ = 2 — 5no- 
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Since we shall consider in this paper only the case of cold quark matter in a magnetic 
background field H, we let T in tend to zero which finally yields 

1 ^ \e-\H °° 

—AV{fi,^,\e,\H) = -J2^-^ dp3E«-^(/^-^^,")(yt^-^^,n) (10) 
This representation will now be used for examining the phase structure of the NJL model. 



3 van Alphen - de Haas oscillations 

Next we will show that in the range of small magnetic fields the expression AV^(/i, S, |ej|i7) 
is an oscillating function of \ei\H which cannot be easily seen from the expression (|TI1|). 
Indeed, equation (|10]) is not quite convenient for extracting relevant physical informations. 
In order to resolve this problem, let us apply Poisson's summation formula 

OO OO 

E a„$(n) = 2 E a„ / cos(27rA;x) dx (11) 

n=0 n=0 "^^ 

to this equation. By integrating first over the momentum variables in (|10[), we find 

3 

1 1 2 OO f\p.\lj\2 

— A\/(/^,S,|e,|if) = - e{ii-Y.)Y.Y.[-^] {S-{a,,k)cosnkh 

-C^{ai,bi) sin Tckbi}, (12) 

where 

S'_(aj, bi) = Si^Tckai) — Sij^kbi), C^{ai, bi) = Cijikai) — Cijikbi). 
The functions C{x) and S{x) are given by Fresnel's integrals 



^ S{x)\ 1 r dy f sm.y\ 



(13) 

\yyx)j \J AT{ J» ^'(/ycosyy 

and = j^q-f^, bi = ^ ' ^ ^ ['-'' Fi^S'lly; after doing some calculations we obtain for 

the oscillating component of AV{fi, S, \ei\H) the expression 

^AVoM, S, \e^\H) = ^^^^7^^ E (^-^] ' {Qinka,) cos{27:ku, + ^) 

TT 

+P{'Kkai) cos{2TikuOi - -)}, (14) 

where uji = fj^^ are the oscillation frequencies. 

Note also that for 1 -C x we have the representations 

P(x) = x'^ -3x-y4: + 
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Equation ( |T4|) is an exact expression for the oscillations of the effective potential of an 
ideal relativistic quark gas at T = 0. The main difference of the relativistic case from the 
nonrelativistic one (discussed in QED by Landau in 1930 ([|l], ||, ^) is the oscillation 
frequency coi = 2\e '\H ^'^ compared with the nonrelativistic expression ujq = ^j-^"^-' . 
Further differences are the dependence of the dynamical quark mass H) on the chemical 
potential /x and the external background field if ([^ 0, |TD|, |TT|). From (|Tip we can easily 
derive the contribution to the magnetization 

l c^(^AKsc(/^,S,|e.|g)) 
(3 dH 

which is obviously an oscillating function of Ui (van Alphen - de Haas effect). 



4 Phase structure of the NJL model 



Next we shall examine the phase structure of the NJL model. For this purpose it is conve- 
nient to rewrite the third component of the effective potential given in (10) in a more useful 
for us form which is completely equivalent to (p!2D. To this end, let us integrate over the 
momentum variable which yields 



■AV{fi,E,\e,\H) 



TC' 



Y,\ei\HY, aJifi - + 2\e^\Hn) 



1=1 



n=0 



X <^ fiJfi^ - S2 - 2\ei\Hn + (S^ + 2\ei\Hn) In ^7 =J 

^Jt,^ + 2\ei\Hn 

It is convenient to consider in some detail the various terms in the gap equation 
1 dV^"\J:,\e,\H) 1 9V^(S) , 1 ^Vo{J:,\e^\H) , 1 dAV{fi,J:,\e,\H) 



+ 



9S 



+ 



9S 



(15) 



(16) 



which leads us to find some new properties concerning the phase structure of the model. 
After some calculations we find for the first contribution 



1 dV^{T.) S 



A", 9S 



F{T) + e{^i-T) 



/iV/i 



(17) 



where FfS) 



A2 



For the second term we have 



1 9K,(S,|e,|iJ) 



Nr. 



In- 



r 



'2ti 



ki ki 1 ki 

H m — 

2 2 2 



:i8) 
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where ki = t^]^, and T(x) is the gamma function. 
Finally, the third contribution is given by 

1 dAV{i2,J:,\ei\H) -2S 



c i=l n=0 



X ^ ^ + In ^ =J }. (19) 
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V-S2-2|e,|i/n) ^S2 + 2|ei|i/n 

It follows from Eqs. (17)-(19) that the gap equation (16) takes the form 

S ■($(/., S,|e,|if)) = S-(0^(S)+0o(S,|e,|/f) + A0(/i,S,|e,|i/)) = O, (20) 

where the function $(/i, S, \ei\H) is represented as a sum of three terms corresponding to 
the R.H.S. of Eq. (16). 

To investigate the phase structure of the model it is convenient to consider in the pa- 



rameter plane (/i, /i), h = ^ max{2\ei\H) , the critical curves Jl{k){h) defined as solutions of 
the equation $(/i, 0, \ei\H) = for the case where the dynamical quark mass S vanishes and 
chiral symmetry will be restored. It is then useful to divide the part of the parameter plane 
in which ^ > 0, h > into regions Qk 

= {(^^ h) : hVk <fi< hVk + 1} (21) 

limited from above by the critical curve fli^k){h). It is obvious that in the region Qq only the 
first term under the summation sign in the contribution from (|19]) is different from zero, in 
Qi the first and second terms are different from zero, and so on. We can show that for fixed 
h the dynamical quark mass S(yu) is monotonically decreasing in /i. In particular, in discrete 
points yUfc = h^/k, k=l,2,..., its derivative turns out to have jumps Ck{h) — Ck-i{h) ^ 
where Ck{h) is defined by 

/dM\ 

C^ih) = -[^] , (22) 

V 9s y ^^^^ 

and for all k one has Ck > Ck+i {k = 0,1,2,...) and Ck < 1. Thus in the region < 
fi < 'Jl(^^^{h) we obtain a nontrivial solution S(/i, \ei\H) ^ 0, whereas for /z > 'Jl(^^){h) the gap 
equation has only the trivial solution. The corresponding behaviour is shown schematically 
in the Fig.l. 



5 Summary 

In this work we have investigated dense cold quark matter in a homogeneous magnetic field 
using an extended NJL-model. Here, we were interested in the case that chiral symmetry is 
broken so that quarks get a dynamical mass E(/x, \ei\H). In particular, we have shown that 
the effective potential of the model is an oscillating function with frequency uoi = frf-m-. In 
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the case \ei\H) ~ fi the relativistic frequency reduces to the well known nonrelativistic 
expression cuq = 0]. Finally, the phase structure of the NJL model in the parameter 

plane (/i, h) and the influence of the oscillations on the /^-dependence of the dynamical quark 
mass \ei\H) have been shortly discussed. 
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Figure 1: Schematical representation of the solution S(/i, \ei\H) of the gap equation (^) 
a function of fi for fixed H with a given critical value /ic = fi(k){h). 
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